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The theory of wave propagation in an inhomogeneous compressible plasma is consid- 
ered. Cylindrical configurations are chosen such that Maxwell's equations, when combined 
with a (single fluid) continuum theory of fluid dynamics, are separable. It is shown that 
genera] expressions for the fields arc superpositions of TM (transverse magnetic), TE (trans- 
verse electric), ;md acoustic waves. For oblique incidence, these three wave types are 
mutually coupled except in special cases. 

1. Introduction 

The scattering of electromagnetic waves in free space from an infinite dielectric cylinder 
is a problem first solved by Lord Rayleigh in 1881 and again by Seitz and Ignatowski in 1905. 
Rayleigh [1918] in a later paper discussed the earlier work and gave some numerical results. 
The complete generalization of the solution to include the case of oblique incidence was at- 
tempted many years later by Wellmann [1937] whose results appear to be incorrect as pointed 
out by van de Hulst [1957]. The complete solution for oblique incidence on a dielectric cylin- 
der was apparently first given only ten years ago [Wait, 1955]. 

An interesting aspect of the obliquity is that incident TE (transverse electric) waves 
will be partially converted to scattered TM (transverse magnetic) waves by the dielec- 
tric cylinder. The converse is also true in the case of incident TM waves. However, if the 
cylinder becomes perfectly conducting, this conversion process no longer takes place. In 
other words, incident TE waves scatter as TE waves and incident TM waves scatter as TM 
waves. The following question is now posed: What happens if the surrounding or ambient 
medium is a plasma instead of free space? Such a situation might occur in scattering of 
electromagnetic waves from metallic cylindrical objects which are immersed in the ionospheric 
plasma. 

It should be self-evident to many that scattering of electromagnetic waves within a cold 
isotropic plasma is not really different from scattering in free space. At least, this is the situa- 
tion when nonlinear terms may be ignored. This follows since the cold plasma, in the absence 
of a d-c magnetic field, may be characterized by a dielectric constant which differs only by a 
complex factor from that of free space. Consequently, an incident TE (or TM) wave should 
still give rise to a scattered TE (or TM) wave provided the cylindrical scatterer is effectively 
a perfect conductor. 

An interesting situation arises when the surrounding medium is a warm plasma. As 
demonstrated below, it is then found that a perfectly conducting cylinder will convert TE to 
TM waves and TM to TE waves. Furthermore, electroacoustic-type modes enter into the 
picture. It is the purpose of the present paper to analyze this question employing rather an 
idealized cylindrical model. Some related questions are also discussed. 
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2. Basic Equations 

The plasma medium is regarded to be a one-component electron fluid. That is, ions are 

neglected in the equations of motion. Also, it is assumed that signals and perturbations are 

sufficiently small for linearized equations to be valid [Oster, 1960; Ginzburg, 1964]. The 

average number density of electrons is n and is taken to be a constant. The pressure deviation 

-> 

of the electrons from the mean is p and their velocity is v. As usual, the electric and magnetic 

-» -» 

fields are denoted E and H, respectively. 

Neglecting any dissipative effects such as collisions, the linearized hydrodynamic equation 
of motion is written 

-> 

& V TV /l \ 

rnn ^-=noeE—vp, (1) 

where e and m are the charge and mass of the electron, respectively. The linearized equation of 
continuity, when combined with the equation of state, is 

u 2 mn V • v= ~-^r (2) 

where u is the velocity of sound in the electron gas. 

Maxwell's equations for the electromagnetic fields in the sourceless plasma are given by 



and 



VX#=-Mo^f (3) 



v/ ft dE { -* , A . 

VXi?==e -ftT+ n o ev > ( 4 ) 



where e and /x are the dielectric constant and magnetic permeability in free space. In what 
follows, all field quantities are taken to vary according to exp (ioot). Furthermore, it is assumed 
that the plasma region under consideration is homogeneous. 

After a certain amount of manipulation of (1) through (4), it is found that the pressure p 
satisfies 

(V 2 +kl)p=Q, (5) 

where 



_co/ cofV 



the 2 

e m 



2 ^^ 
and ooq^ 



On the other hand, by making use only of (1) and (4), it is found that the magnetic field is 
obtained from 

H=-VXA, (6) 

where A satisfies 

( V *+k 2 e )A=0, (7) 

and where kl^e^oy 1 and (e/e Q ) = l — (co /w) 2 . As seen by its definition, e is the dielectric constant 
of the plasma. In the case of a cold plasma, the electron acoustic velocity u tends to zero so 
that the acoustic wave number k p becomes infinitely large. 

It is of interest to note that the pressure deviation p is related to the static pressure p and 
density perturbation n by the adiabatic relation 

p=yp nln , (8) 
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where y } the classical ratio of specific heats, takes the value 3 for 1 degree of freedom in the 
gas. Thus, in terms of the kinetic temperature T, it is found that 

P=3kTu, (9) 

where k is Boltzmann's constant. In terms of the root-mean square electron velocity u, 

p=mu 2 n. (10) 

For the remainder of this paper, we shall refer only to the pressure deviation p but one should 
keep in mind that this is related directly to the density perturbation of the plasma [Cohen, 
1961]. 

On combining (1) and (4) for the harmonic time dependence it is easy to see that 



and 



E=J- (vxfl)- ,f 2 , Vp (11) 



(VXff)-^— Vp. (12) 

iumn e 



3. General Field Representations 

In the present paper, we are particularly interested in obtaining representations which 

are applicable to cylindrical coordinates (p, </>, 2). Thus, A and p are first expressed in a Fourier 
representation in the manner [Wait, 1959] 

4=£ f a n e- ihz dhe- in + (13) 

n = -oo J -co 

and 

_f_ co /»+co 

P= 2 I3 fi e- ihz dhe- in +, (14) 

where a n and p n are the respective transforms. 

In operational notation, equations (13) and (14) may be written 

A=Yt n (13') 

and 

P=rp n (140 

where r signifies multiplication by exp [— ihz— irut>] and subsequent integration over h, and 
summation over integer values of n. 

In view of (5) and (7), it is seen that a n and ft n satisfy 

(V?+^ e )a n =0 (15) 

and 

(Tt+OP»=0, (16) 

where V 2 t is the two-dimensional (i.e., transverse) Laplacian operator, while 
and 

1/p r== /^p lb . 
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In cylindrical coordinates (p, </>, 2), it is well known that 

' pd P p dp^p'dcjy 2 ~pdp p c)p p 2 ' 

Therefore, solutions of (15) are cylindrical Bessel functions Z n (u 6 p) of order n and argument 
u e p. Similarly, solutions of (16) are cylindrical Bessel functions Z n (u p p). In the subsequent 
development, it proves to be convenient to use the Bessel functions J n and the Hankel function 

H^ as the two independent solutions. Furthermore, rather than working with the vector 

-> 
potential A, it is more convenient to employ the two scalar quantities E 2 and H 2 in a manner 
similar to that used for dielectric media [Wait, 1959]. 
With some consideration, we now find that 

E z =T[u 2 ^ n H n +ulA n J n ]+-^ r T[hc n Hn+hC n J n ], (17) 



and 
where 



H z =T[u.b n H n +u.B n J n ] 9 (18) 

p=T[c n H n +C n J n ], (19) 

H n =m\u e p),H n =m\u v p), 

A 

Jn = Jn(u e p), J n = Jn(u p p), 

and a n , A n , b n , B n , c n , C n are coefficients which do not depend on the coordinates. 

Using the basic equations (1) to (4), it is a straightforward matter to obtain, from (17), 
(18), and (19), general expressions for the transverse components of the fields. Thus, we find 

E,= T [iwh ^?»--5* a n H n +iwB H ^- n Jt An j 1+-^- T[nc n H n +nC n J n ], (20) 

L Op p Op p J moo ep 

(21) 

H* = T |~-- JA-^ a n ^=-2* B n J n -^- A„ ^], (22) 

L P Mo^ Op p /x co Op J 

H P =T \-ihb n ^+^L a n H n -ihB n ^+^1_ ^j "]. (23 ) 

L Op /x cop Op /x cop J 

The components of the velocity v are now found easily by using (12). Thus, for example, the 
radial component is given by 



Vp= 2" 

orme 



/l dH z dff,\ dp . 

\P d(j> dz ) ioomrioe dp 

or, explicitly, 

r [- ( - VXH-^ klB n J n +^ a n ^+^ An ^1 

LP P Mo^ Op Mow Op J 
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4. Treatment of Cylindrically Stratified Media 

The representations developed above are particularly suitable when dealing with cylin- 
drically stratified media. For example, we could deal with a column of a compressible plasma 
whose properties varied only in the p direction. The idea would be to represent the region by 
a finite number of concentric cylindrical homogeneous regions. In a typical region bounded 
in p by p t and p i+1 , the elect roacoustic velocity is u t and the electron plasma frequency is (co )*. 
In the region p m to p i+2 , the corresponding plasma properties are then designated n ni and 
(o> )h-i. In each region, the field representations would have the form given by (17) to (25). 
Thus, in the iih region, the six unknown coefficients are a in , A in) b im B in , c in , and ( \ n . In t lie 
(^+l)th region, there is a corresponding set of six coefficients. In principle, the boundary 
value problem is solved if six linear equations can be obtained from the boundary conditions 
at p=pu 

Under the rather idealized assumption that the waves do not deform the boundary, 
kinematic considerations require the continuity of the normal fluid velocity. Thus, 

(v P )i=(v P )i+i atp=p z . (26) 

It is clear that, if this condition were not satisfied, the fluids would separate at the boundary 
as pointed out by Field [1956]. 

A further hydroclynamic-type boundary condition is obtained requiring that the pressure 
is continuous. Thus, as indicated by Field [1956], 

Pi^Pi+i atp= Pi , (27) 

which is a requirement imposed by force equilibrium at the boundary. The idealized kinematic 
and hydrodynamic boundary conditions have also been u^ed by Yildiz [1963] and Johler [1964] 
who considered scattering from a compressible plasma sphere. 

Four additional boundary conditions are readily obtained by imposing the usual electro- 
dynamic property that tangential electric and magnetic fields are continuous. Thus, 

(E z ) t =(E 2 ) i+1 ' 

(E<t>) i = (E<t>) i+i 
(H z ) t = (H z ) i+ i 

(H (j) ) i =(H 4> ) i+ i 



at p= Pi . (28) 



For P+l homogeneous regions, there are only P interfaces. Thus, it would appear, in 
general, that insufficient boundary equations are available. Fortunately, however, additional 
physical considerations allow us to dispense with some of the coefficients in the outermost and 
the innermost cylindrical regions. For example, if the source is in the outermost (Oth) homo- 
geneous region, the coefficients (A n ) , (B n ) 0) and ((7 n ) would be known since they are specified 
by the form of the incident field. The scattered field in the outermost region would be described 
entirely by the coefficients (a n ) 0} (6 n ) , and (c n ) . On the other hand, in the core or innermost 
(P+lth) homogeneous region, the coefficients (A n ) P+lJ (B n ) P+1 , and (C n ) P+1 are zero since 
the fields must be finite at p=0. 

5. Application to Rigid Perfectly Conducting Cylinder 

As an illustration of this method for handling cylindrical boundary value problems in 
compressible plasma, we shall consider a highly idealized situation. The cylindrical boundary 
at p=a is regarded as perfectly conducting and rigid. The exterior region p>& is taken to be 
a homogeneous compressible plasma with properties u and co . With these simplifying and 
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rather restrictive assumptions, it follows that the required boundary conditions are 



v p =0 
E z =0 



at p=a. 



(29) 



The same conditions have been employed by Seshadri et al. [1964] in the study of scattering 
of normally incident plane waves from a perfectly conducting cylinder immersed in the com- 
pressible plasma. Using (17), (20), and (25), it now readily follows from (29) that, in matrix 
form, 



TL 



u 2 e H n 

nh 
a 

k 2 e hu c 

/X CO 



H' n 





ijJL O>Hn 

in g k 2 



Hn 



i g hH n 

A 

ing H n 

€qU p 



A / 

ioomn e 



b n 



-u 2 e J n A n +0—i g hJnCn 
nh 
a 

k 2 e hu ( 



J n A n —in o)JnB n —im g J n C n 



fl O) 






iumn e 



J'nC n 



(30) 



where the arguments of the Bessel functions are evaluated at p=a and where g = «-• 



The 



column matrix on the right of the previous equation contains the driving terms which are 
specified by the form of the incident field. Thus, in principle, the boundary-value problem 
is solved since a n , b n , and c n are now expressed in terms of known quantities. 

A special case of some interest is the cold plasma. The transition to this case is obtained 
by letting the velocity u-^0, whence k p -^co . From (30), we readily find that 



and 



_ Jn A 

TJ *±n 
tin 



/) J n r> 

Un jjr J-^n~ 

t± n 



J n (u e a) 



H™{u e a) 



Jn(u e a) 



A n 



B tl 



(31) 



(32) 



H™'(u e a) 

The first of these equations indicates that if the incident wave is purely TM and characterized by 

ET=T[A n J n ], HT=0, (33) 

then the secondary field is also purely TM and given by 

E z =r[a n H n ], flj=0, (34) 

E } =T [-* a n H n ]> H° = T [-^ a, <f } (35) 



Ei^i^^} m-r[^aA] 



(36) 



(where the arguments of the Bessel functions are ua). 
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Similarly, (32) indicates that if the incident wave is purely TE and characterized by 

Hr=T[B„J n ], E™=0, (37) 

then the secondary fields are also purely TE and given by 

Hl=T[b n H n ], E s 2 =0, (38) 

//;=r [-^ b n Hn~\> e;=t [i^b n ^], (39) 

H;=T \-ihb n ^], El=T [-» b n H^ (40) 

As seen above, scattering of electromagnetic waves from a perfectly conducting cylinder 
in a cold plasma is identical in form to scattering of electromagnetic waves from a perfectly 
conducting cylinder in free space (e.g., as discussed in Wait, 1959). The only change required 
is the replacement of the dielectric constant e of the plasma with that for free space € . 

Another interesting special case occurs when we let the charge density in the medium ap- 
proach zero. Thus, setting e=0, it is easily found from (30) that 



A 



p __ Y_n_ p __ «/ n\ U p a ) p (A\\ 

where u p = (k 2 p —h 2 ) l/2 and k p =(a/u in terms of the acoustic velocity u for an uncharged medium 
of n particles, of mass m, per unit volume. For this case, an incident acoustic wave, whose 
pressure perturbation is given by 

r c =T[C n J n (u pP )} } (42) 

leads to a scattered acoustic wave of the form 

r c =T[c n m\u p p)}. (43) 

This is the well-known result for pure acoustic scattering from a rigid cylindrical obstacle 
[Morse and Feshbach, 1953]. 

We see from the above that our general solution contains both pure electromagnetic and 
pure acoustic scattering as special cases. In the general case, it is easy to see from (30) that a 
pure incident TM electromagnetic wave (i.e., A n ^0, B n =C n =0) will generate scattered TM 
and TE electromagnetic waves in addition to a scattered acoustic wave. In a similar fashion, 
a pure incident TE electromagnetic wave (i.e., 5„^0, A n =C n =0) will scatter TM and TE 
electromagnetic waves and the acoustic waves. Of particular interest is the fact that a pure 
incident acoustic wave (i.e., C n ^0, A n =B n =0) will generate both the TM and TE electro- 
magnetic waves in addition to the acoustic waves. 

Consistently, in the development given here, the field functions are preceded by the 
Fourier-summation operator r. Using this artifice allows us to consider quite general types 
of incident fields. For example, if an obliquely incident plane TE wave is considered, 

E'z c =E sin exp [ik e p sin 6 cos <j>—ik e z cos 0], (44) 

where 6 is the angle which the wave normal subtends with the z axis. Employing a well- 
known addition theorem [McLachlan, 1934], it is found that an alternate representation is 

E™=E sin X) i n Jn(k e sin 6) exp [—in<t>—ik e z cos 0], (45) 

n = — oo 
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which may be written in the form 

ET=T[u 2 e A n J n (u eP )] (46) 

if u%A n =E (sin 0)i n , h=k e cos 6, and u e =k e sm 6. The operator for this case of plane wave 
incidence implies multiplication by exp [—incj)—ik e z cos 6] and subsequent summation over 
integer values of n from — oo to + °° • 

The total fields in the compressible plasma exterior to the rigid perfectly conducting cylin- 
drical object are thus given by (17) to (25) with the above meaning to the operator r. The 
coefficients a n , b n , and c n are given by (30) where B n =C n =0 and A n is defined as above. 

The explicit form of the fields for obliquely incident plane TM electromagnetic waves or 
a plane acoustic wave is obtained in a manner almost identical to the above. The reader 
should have no difficulty satisfying himself that this is the case. 

A case of some practical importance is when the rigid cylinder considered above is excited 
by an electric dipole located at (p , </> , z ) where p ^>a. The electric dipole of length ds with 
current / is taken to be oriented in the ^-direction (i.e., parallel to the cylinder axis). From 
earlier work, we know that [Wait, 1959] 

ET=i^- £ € -<»<*-*» (HZ\u Po )J n (up)exp[-ih(z-z )]dh, (47) 

O7reco n = — oo J 

— 00 

provided p<^po- (For p^>p , we merely interchange p and p in the above representation.) 
Thus, we may write (for p< p ) 

ET=T[A n J n (u eP )}, (48) 

where 

Tel* 

A n =i tP 1 H% ) (u e p Q )e in +°e ih2 ° (49) 

and r signifies multiplication by exp (—ihz—in(t>) and subsequent summation over n and 
integration over h. 

The total fields are now again obtained from (17) to (25) with this general meaning at- 
tached to T. Thus, the fields are represented as a Fourier integral and an infinite sum. The 
formulas can be greatly simplified in the far field where both k e p and t p p>>l. Then the 
scattered field integrals may be asymptotically approximated in the manner [Wait, 1959] 

-| fF(h)H^(up) exp [-ih(z-z )]dhg±exp (imr/2)F(k sin 0) exp (-ikR)/R, (50) 

— CO 

where u=^ e or u p and k=k e or k P , which is valid when the function F(h) is relatively slowly 
varying. Using this asymptotic relation, formulas for the radiation patterns of the dipole in 
the presence of the cylinder are readily obtained. 

When the observer is near the cylinder, it is necessary to consider the influence of the pole 
singularities of the integrand. For example, F(h) may have a finite number of complex poles. 
Some of these may lead to residue contributions which are identified with axial surface waves. 
This would be an important subject for further investigation. It will not be pursued here. 
The interested reader may consult a fundamental paper by Samaddar and Yildiz [1964] who 
consider the excitation of axially symmetric (i.e., n=0) surface waves in a compressible plasma 
cylinder surrounded by free space. 

6. Two-Dimensional Situations 

Another important special case is when the fields do not depend on z [Wait, 1964a]. This 
occurs, for example, if the excitation is a normally incident plane wave or a cylindrical wave 
emanating from a uniform line source. In the general formulas listed from (17) to (25), this 
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means that h is set equal to zero and the operator r now is defined by 

P=T(3 n = Tj P n e- in +. 



(51) 



and 



It is evident that for z independent fields, (30) be partly decomposed to the set 

H n a n =-J n A n (52) 



ifiQ uk e IT n in <P II n 



in£Jc e 



e okp jji 



loomrioe 



-i/jL <dk e J' n B n —iii <P J n C n 



i(jomn e 



(53) 



Here, the arguments of the Bessel functions J n and H n are k e a while the arguments of J n and 

A 

H n are k p a. 

It is evident from the above that a normally incident TM wave will scatter only a TM 
wave even though the plasma is compressible. However, a normally incident TE wave will 
scatter both a TE wave and an acoustic-type wave. To illustrate this point more clearly, a 
specific case is considered. The incident TE wave is defined by 



Hr=H exp (ik e p cos 6) = II 2 i n e- in +J n (k e p),P l] 



--0. 



(54) 



By utilizing the matrix equation (53), it is found that the secondary TE electromagnetic wave 
is obtained from 

Hf=H S i n A T M 2) (k e p)e- in *, 



where 



K~- 



(k e a) (fcX^nfea)#f'^ 



(k e a) (k p a) HT ' (k e a) H? ' (k p a) -n 2 (u /a,) 2 H™ (k e a) H? (k p a) 
It also follows from (53) that the scattered acoustic wave is to be obtained from 

p=p sc =II 2 i n PnII { n ) (k P p)e- in \ 

n = — a 

(2in/if) (o>o/") 2 (jnu/e) 



where 



Pn = 



{k e a) (k p a) II™ ' (Jc e a) H* ' (k,a) - n 2 (coo/^H™ (k e a) H? (Jfcyi) 



(55) 
(56) 

(57) 
(58) 



In writing (58) in the above form, the well-known Wronskian relation for Bessel functions has 
been utilized. Equations (55) and (57) agree with the results derived directly by Seshadri 
et al. [1964]. Their analysis was restricted at the outset to normal (i.e., perpendicular) 
incidence. 

It is easy to see that (56) may be rewritten in the form 



where 



and 



A n =A c n [l— A n ], 

d n \K e d) 



A£= 



H™'{k e a) 

n 2 (W") 2 JnM H?M 



1-A. 



(k p a) 2 (k e a) 2 J' n (k e a) H^'(k p a) 
n 2 (W«) 2 m\k e a) H™{k p a) 



(59a) 
(59b) 

(59c) 



(k p a) 2 (k e a) 2 H™\k e a) H™'{k P a) 
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This shows that for a cold plasma (i.e., k p a-> co), A w — A£, which is the appropriate value for 
scattering (at normal incidence) of a plane wave from the perfectly conducting cylinder im- 
mersed in a cold plasma whose wave number is k e . If either u/c or w /a> is sufficiently small, 
the correction factor A n itself may be very small. To within this approximation, it is not 
difficult to show that 

A _ (2n 2 /T)(a>o/o>) 2 (u/c) 



(k e ayj' n (k e a)H™'(k e a) 



(60) 



where, as usual, k e =k(\— cog/co 2 )' in terms of the free space wave number k. In most situations, 
the factor A n is exceptionally small. For example, in the case of ionospheric plasma u/c is of 
the order of 10~ 3 . 

7. Concluding Remarks 

On the basis of the analysis in this paper, it is evident that TE, TM, and acoustic waves 
are coupled in a cylindrically stratified compressible plasma. However, when the fields are 
independent of the axial coordinate, the TM waves become decoupled, although the TE and 
the acoustic waves are still coupled. Also, in the limiting case of a cold (i.e., incompressible) 
plasma, the acoustic wave is not present but there may still be coupling between the TE and 
TM electromagnetic waves. 

In the analysis in the present paper, the presence of a superimposed d-c magnetic field is 
ignored. As indicated by Field [1956] and much more recently by the author [Wait, 1964b], 
this will cause additional coupling. Furthermore, the influence of nonlinearities in the medium 
will greatly complicate matters [Ginzburg, 1964]. 
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